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We demonstrate how Bose-Einstein correlations emerge from the correlations of 
fluctuations allowing for their extremely simple and fast numerical modelling. Both 
the advantages and limitations of this new method of implementation of BEC in 
the numerical modelling of high energy multiparticle processes are outlined and 
discussed. First applications to description of e"'"e~ data are given. 

1 Introduction 

Bose-Einstein correlations (BEC) between identical bosons are since long time 
recognized as indispensable tool in searching for dynamics of multiparticle pro- 
duction processes because of their potential ability to provide some space-time 
information about them ^ However, serious investigation of such processes 
can be performed only by means of involved numerical modelling using for this 

Surpose one of the specially designed Monte Carlo event generators (MCEG) 
. Their a priori probabilistic structure prevents the genuine BEC to occur 
because they are of purely quantum statistical origin. The _biea|-one can do 
is either to change accordingly the output of such MCEG's oB-eHZI to provide 
for the necessary bunching of identical bosonic particles (mostly tt's) in the 
phase space or to try to incorporate somehow their bosonic character (or at 
least some specific features distinguishing the bosonic and purely classical 
particles) into the MCEG itself, i.e., to construct the MCEG providing on 
the output particles already showing the bunchiijLE_mentioned above l3. Our 
proposition discussed here (based on our works EHij) can be in this respect 
regarded as belonging to the first category but, at the same time, it uses alsD 
ideas explored in □ and, to some extent, can be regarded as improvement of □ 



"The "infinite" number of references on BEG-«vill be projected here on some selected number 
only, starting from the pedagogical review a and followed by a selected number of more 
recent, mostly complementary in their approach to BEC, reviews p. 
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(at least in what concerns the numerical performance). 

Whatever one is doing the final objective is always the same: to reproduce 
the characteristic signals of BEC obtained experimentally, which for the case 
of 2-particlc BEC means that the following two-particle correlation function, 



C2{Q^\P:-P,\) 



(1) 



defined as ratio of the two-particle distributions to the product of single- 
particle distributions increases towards C2 = 2 when Q approaches zero. No- 
tice that ([^) depends explicitly on the measwred momenta {pi,pj) of the (like) 
particles, the searched for space-time information can be obtained only when 
treating it as a specific Fourier-transform of the distributions p{r) of the pro- 
duction points, in which case C2 is usually (schematically) written asEllj 



C2iQ) 



dx p(r) 



ip(o)r 



(2) 



It allows (in principle) to translate the (observed) width of the peak in C2iQ) 
into the (deduced) size of the region of emission p{r) ^ From that point of 
view the first approach to the numerical modelling of BEC mentioned above 
is addressing (in one or other form) directly the (^ whereas second (and our 
as well) is providing only (^ which can be later analysed (much in the same 
way as the experimental data are) to obtain the information on p(r) 



2 BEC - quantum-statistical approach 

In what follows we shall treat the BEC as arising because of correlations 
of some specific fluctuations present in physical system under consideration 
(known as photon bunching effect in quantum optics c3 where similar corre- 
lations are also known under the name of HBT effect). Notice thatJjjic main 
ingredient of C2 is the correlator (nin2), which can be written as EHIj || 

{nin2) = {ni){n2) + ((ni - (ni)) (^2 - ("-2))) = {ni){n2) + pa{ni)a{n2). (3) 



''In reality all this is, of course, much more complicated El'0 but this will not be our point of 
interest here. Notice only that what is being deduced in this way is not the whole interaction 
region but rather the region where the like-particles with similar momenta are produced. 
We could call it the elementary emitting cell, introducing in this way the idea used later 
on in this talk. 

'^Here cr{n) is dispersion of the multiplicity distribution P{n) of produced secondaries and 
p is the correlation coefficient depending on the type of particles produced: p = +1, —1, 
for bosons, fermions and Boltzmann statistics, respectively. 
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Therefore the two-particle correlation function (|^) is entirely given in terms 
of the covariances ^ explicitly showing its stochastic character: 

In eq.(^ above C2{Q) is just a measure of correlation of fluctuations present 
in the system under consideration, which is maximal (i.e., C2 = 2) whenever 
fluctuations are maximal (i.e., for <T{n) = (n), what happens for the geometri- 
cal distribution of the produced particles, which in turn happens in a natural 
way where they are bosons and as such show maximal tendency to bunch 
themselves in the same state (or cell) in the phase space). 

This feature is the cornerstone of the only MCEG providing particles al- 
ready showing C2{Q) > 1 without any additional procedures]^. Those who 
would like to start with the synuactrization of the corresponding multiparti- 
cle wave function should realize 113 that such symmetrization is equivalent to 
the change from the Maxwell-Boltzmann statistics (classical, distinguishable 
particles) to the Bose- Einstein one (indistinguishable quantum particles). To 
do this one has to select as independent subsystems not the individual parti- 
cles but rather the groups of particles in the consecutive states. In effect the 
original Poissonian distribution of particles goes over to the geometrical one 
mentioned above. This can be visualised in a best Wjay on the following simple 
example, which will also represent the main ideas of l3|^. Suppose that mass M 
(in rest) is going to hadronize (for simplicity for one kind of neutral particles 
with mass m each). When one is selecting each time a particle with energy 
Ei according to the simple statistical distribution P{Ei) — exp {~Ei/T) until 
the whole M is used up, one gets multiplicity distribution of the produced 
secondaries, P(n), in poissonian form (with < n > depending on the parame- 
ter T). Suppose now that one is changing algorithm in a following way: after 
selection the first energy, Ei, one adds with a probability P to the particle 
chosen in this way another particles of the same energy Ei and does it until 
the first failure (i.e., until the random number selected is greater than P). 
Then one is selecting a new particle with energy E2 (i.e., in fact one is select- 
ing a new energy E2) and repeats the above procedure again and again, until 
the whole mass M is used up. It is straightforward to realize that what one is 
getting is the number of poisson-like distributed cells (with {riceiis) depending 
on T containing each a number of geometrically distributed secondaries (with 
(n) in a given cell given by parameter P). Taken together this convolution 

''The nearest in spirit to is the unpublishnd attempt presented in0. The two others are 
providing BEC in not so natural way as Q. n 
Those interested in details should consult Q directly. 
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of poissonian and geometrical distributions results in well defined Negative 
Binomial (NB) distribution of the produced secondaries ^. However, in or- 
der to get C2{Q) with a characteristic width AQ ^ 1/R (corresponding to 
a "radius" R) one has to allow for a spreading of energies i?^*'' of particles 
belonging to the fc-th cell, A^;. Let us stress at this point that it is precisely 
this spreading which translates finally into the dimensional parameter R being 
the main subject of interest when interpreting experimental results for C2. 

The MCEG proposed in U has the same structure with particles selected 
according to Bose-Einstein distribution, P{Ei) ~ exp [nj — i^^) /T] (n^ is 
their multiplicity and Ei are their energies, the two parameters, "temperature" 
T and "chemical potential /i correspond to the previous T and P) and the 
like-particles are located in the cells of a fixed size 5y (in the longitudinal, 
i.e., rapidity y space), which is the third parameter corresponding to the A^; 
above). 

It turns out that in this approach one gets at the same time both the 
correct BEG pattern (i.e., correlations) and fluctuations (as characterized by 
the observed intermittency pattern) This is very strong advantage of this 
model, which is so far the only example of hadronization model, in which Bose- 
Einstein statistics is not only included from the very beginning on a single 
event level, but it iS|-als£j-.properly used in getting the final secondaries. In 
all other approaches □laQQ at least one of the above elements is missing. The 
serious shortcoming of method u are, however, numerical difhculties to keep 
the energy-momentum conservation as exact as possible and its limitation to 
the specific event generator only^ 

3 BEC - our proposal: general ideas 

In we have proposed a new method of numerical modelling of BEC, which 
we shall -shortly introduce now (some of its elements were already formu- 
lated in □). It was thought as additional element of any MCEG describing 
multiparticle production processes. So far it was tested only_pn some spe- 
cial hadronization scheme proposed by us (the CAS model Ea) and on the 



■''This can shed a new hgbjt on the NB distributions and its different applications as discussed 
during this conference |^ with, for example, our cells being the equivalent to clans in the 
standard approach to iWB distributions. Similar concept of elementary emitting cells has 
been also proposed in Q 

f Notice that in our simple example, which apparently did not show these shortcomings, we 
were considered only one type of produced secondaries. When one attempts to incorporate 
all charges and to keep the energy-momentum conserved, our example starts to be also 
intractable in practice. 
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JETSET model for e+e~ processes (to be shown below). We aimed at fast 
and universal algorithm providing BEC among the secondaries produced by 
a given MCEG already at the event-by-event basis. 



Our way of reasoning was as follows: the BEC is a quantum mechanical 
phenomenon whereas all MCEG are of classical character. Therefore to mimic 
BEC one has to resign from a part of information provided by the MCEG, 
which usually gives us energy- momenta, spactj^time positions and charges (Qi) 
of the produced secondaries. All methods BHq'EJQ'I] are in one or other way 
changing the first twc|^ We propose to keep them intact and to change, 
instead, the originaLcharge assignment provided by the MCEG. This is done 
in the following wayEiJ. Suppose that our MC event generator provides us with 
N{+), N{—) and A^(0) of positive, negative and neutral particles, uniformly 
distributed and showing no BEC pattern (cf. Fig. 1, left panel). We change 
now their charge allocation (keeping the same N{+), N{~) and N{0)) and 
getting the picture shown in the right panel of Fig. 1. The like charges 
are in visible (albeit strongly exaggerated) way bunched (correlated) together 
leading to signal of BEC. What we have done is the following: (a) we have 
resigned from the (not directly measurable) part of the information provided 
by event generator concerning the charge allocation to produced particles 
and (b) we have allocated charges anew in such a way as to keep the like 
charges as near in phase space as possible (keeping also the total charge of 
any kind the same as the original one). In this way we have formed abjects 
which we shall call in what follows elementary emitting cells (EEC) each 
containing only particles of the same sign and belonging to the same state. In 
Fig. 1 they are shown as separate in the phase space but in principle they can 
overlap. Particles belonging to such a cell are supposed to be in the same state 
(the fact that their momenta differ from each other reflects only the natural 
spreading of such state in the momentum space and precisely this spreading 
is the source of the analogous complementary spreading in the position space 
and eventually in the observed structure of the C2{Q))- 

The actual implementation of our proposition (algorithm) is presented 

''This is especially visible lig works claiming to describe BEC from the very beginning in a 
quantum mechanical way 
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in detail in ll^l. Here we shall discuss some further points of the proposed 
approach and present new results of its application to the known algorithms 
for hadronization applied to the e+e^ data. We would like to stress the 
point that the basic entity here is not so much the hadronizing source but 
the EEC mentioned above and C2{Q) are in our approach directly sensitive 
to the number of such cells and to their mean occupation. There is no direct 
dependence on the "size of the hadronizing source" |. On the other hand our 
approach contains automatically BEC of all orders (in practice, the highest 
order accounted for is given by the highest multiplicity in a given EEC one 
can reach). 



4 BEC - our proposal: numerical results 

To illustrate how our prescription works we present in Fig. 2 sample of results 
using as initial source oLsecondaries the CAS hadronization mod el developed 
by us some time ago E3 and the standard JETSET model EJ. This very 
preliminary and limited in scope attempt (for example, for simplicity only 
direct pions were produced both by the CAS and JETSET MCEG's) shows 
in a clear way the main points of our method mentioned before. The EEC's 
were formed according to the algorithm given in 113 using as parameter the 
probability that a new selected particle will joint the particular EEC which 
is being actually filled. This probability was kept constant in the case of 
CAS (and equal P, as indicated in Fig. 2(a)) whereas in the case of using 
JETSET it was given by P = exp(— S/T) (where E is the energy of the 
particle considered at given moment) with value of parameter T as indicatec 
in Fig. 2b. In the case of CAS the, so called, "two split- type sources" (cf. 
for details) were used[]. Together with fits to the sample of DELPHI data 
are shown the corresponding distributions of EEC's, P{Nceu), and particles 
in such cell, Pijipart)- 



'Notice, that the "size" -R (or its variants) discussed in all known analyses of BEC (both 
data and formulas) are representing at first the distance between the two emitting points, 
R = n — r2, and therefore are only indirectly sensitive to the "true" «ize of the source. 
Nevertheless, customary they are referred to as the "size of the source" El. 
^ It means that original source were supposed to consists of two sub-sources of equal mass 
each, which were cascading independently, and our algorithm was then applied to all parti- 
cles without discriminating which source they were coming from. In effect one is obtaining 
much more dense source (i.e., partjctes are packed more closely in the phase space) than in 
the oftse of the single cascade only 1121. The slope of C2{Q) turns out to be very sensitive to 
this tj. The opposite situation, in which particles from different sources are supposed not 
to show BEC, is referred to as "indep-type spwpces" and can, for example, easily explain 
the BEC puzzle for W^W~ production data [j. 
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gure 2. Application of our method to fit DELPHI datpftj on e+e~ annihilation using CAS 
t3 (a) and JETSET7.4.10 (with standard parameters) Ej (b) hadronization models. Panels 
(c) and (d) show the corresponding distributions of EEC, P{N^^ii), and distribution of the 
like-charge particles allocated to such cell, P{npart)- The best poissonian and geometrical 
fits to the, respectively, PiN^^n) and P(npart) are also shown. 
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Notice that, with a very good accuracy, distribution of particles in a sin- 
gle EEC is of geometrical type (i.e., according to our discussion before, corre- 
sponding to boapnic particles). This follows directly from the construction of 
our algorithm Ii3 (one is adding new particles to the EEC with probability P 
until the first failure, in such a case in the ideal situation (upart) = P/ i^ — P))- 
The strength of BEC will in our case depend on how much (on average) 
particles are allocated to EEC and how many EEC's one has (on average) 
in an event. One can see in Fig. 2 that the first number is rather low, 
{npart) ^ 1-2 1.3, the probability to get 3 or more particles in a single EEC 
is of the order of 1% only (and less). It means that multi-particle BEC are 
not very prominent feature under normal circumstances and, if at all, should 
be probably looked for in a ijeally very high multiplicity events (which, on 
the other hand, are very rare Ell). The second number is of the order of ~ 7. 
Notice that the EEC's are distributed essentially in a poisson-like manner. It 
is interesting then to notice that, according to what we have already men- 
ned before, it means that the total multiplicity distribution is of NB type 
Actually, this is the distribiition provided hji_the original MCEG we are 
working with (in our case CAS El and JETSET El), which is then reproduced 
by the action of our generator modelling the BEC (as it should be)^. 

5 Summary 

Our presentation of BEC is surely not the orthodox one, i.e., we do not see 
BEC as the immediate source of the information on the space-time character- 
istic of the hadronizing object. We approach the problem from the quantum 
statistical point of view, stressing more the behaviour of the like-charged 
bosons in the momentum-space (which, after all, is what is really measured!) 
necessary to be incorporated in the MCEG whenever one aims to model the 
true BEC. Our particles feel the BEC only when they are in the same state 
(it allows for the sizeable differences in their momenta because of the spread- 
ing of momenta of the wave packets representing such states). It means then 
that our primary objects are EEC's rather, than the whole emitting hadronic 
source itself. And it is the property of the average EEC, which is responsible 
for the finally observed structure of the C2{Q) function. 

Most of our presentation deals with a kind of universal algorithm 1121, which 

'"It is worth to stress at this point that, according to our discussion in|^, a constant number 
k of EEC with geometrical distribution of particles in each of them leads to NB (and in the 
limit of large k to Poisson) distribution whereas binomially distributed EEC's, with k limited 
for some reaperi, leads to the so called modified negative birtemial (MNB) multiplicity 
distributions E3 characterised by oscillating cumulant moments 
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should be suitable as a simple addition to (almost) any MCEG |]. However, it 
is very likely that such approach will fail when subjected to the more serious 
scrutiny than it is done so far (in the sense that its action will have to be 
connected with a re-parameterization of the original MCEG in order to get 
the right resul±s_i.e., to correct for deviations introduced by the BEG imple- 
mentation Bu'EIeIi]). In such case, the only solution left seems to be to start to 
construct MGEG from assuring that it models in a correct way the bosonic 
character of producecLsecondaries, i.e., to follow and improve the program 
started with the work □ along the lines presented at the beginning of this talk. 
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